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ABSTRACT 


This report presents the results of a validation study of data obtained 
by the Stratospheric Aerosol and Gas Experiment II satellite experiment (SAGE 
II). Preliminary SAGE II data have been available for the period October, 
1984 to May, 1985. In addition, the results of two correlative experimental 
measurement series have been studied in detail, as well as climatological 
data obtained by other techniques, including ground-based and airborne 
lidar. The study shows the SAGE II data to be of great potential value to 
studies of the microphysics of stratospheric aerosols, the chemistry of trace 
gases and stratospheric dynamics. A small number of unidentified errors in 
the current preliminary data set are described. These will be removed from 
the next version of the data set which is anticipated to be of archival 
qual ity . 
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1. INTRODUCTION 


1 . 1 Background 

The Stratospheric Aerosol and Gas Experiment II (SAGE II) was launched 
aboard the Earth Radiation Budget Satellite (ERBS) in October 1984 (McMaster, 
1986). It is the fourth in a series of satellite-borne stratospheric aerosol 
monitoring instruments flown by NASA. The SAGE II instrument is a 
seven-channel solar radiometer which measures solar ocultation at the chosen 
wavelengths for 30 events on each day. The specific objectives of the SAGE 
II experiment are to: 

1. map vertical profiles of stratospheric aerosols, ozone, nitrogen 
dioxide and water vapor with a 1 km vertical resolution and determine high 
altitude cloud coverage from 80°N to 80°S; 

2. study the seasonal and global variations in these constituents so as 
to improve our understanding of the Earth's radiation budget and climate; 

3. utilize the measurements to study global circulation, transient 
stratospheric phenomena (such as volcanic inject ion) and sources and sinks of 
the individual species; and 

4. utilize the measurements to study the atmospheric chemistry of the 
gaseous species and the microphysics relating to the aerosols. 

Before data from SAGE II or other similar experiments may be used for 
the scientific objectives listed above, it is necessary that they undergo a 
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validation process. During this process, the data are compared with 
preexisting models of stratospheric aerosols and gaseous constituents, as 
well as compared with correlative data sets obtained by other means. The 
statistical nature of the data product must be examined both for internal 
consistency and for external agreement. This report describes work carried 
out on the available SAGE II and related data sets as part of this validation 
process . 

1.2 SAGE II and Associated Data Sets 

The basic SAGE II measurement is of solar radiance in the seven channels 
listed in Table 1.1 (McMaster, 1986; Mauldin et al,, 1985) made during the 
spacecraft sunrise and sunset periods. The location of the measurements is 
determined by the spacecraft orbit. This provides 15 sunrise and 15 sunset 
events on each day. On a specific day, the sunrise events occur at 

approximately the same latitude as do the sunset events, the latitude of the 
sunrise events normally being different from that of the sunset events. 
Successive sunrise or sunset events are separated by approximately 24^ of 
longitude. The event latitudes cycle from one extreme of latitude to 

another, a cycle covering about 120^ of latitude and taking about 5 weeks to 
complete. The extremes of latitude vary somewhat with season but lie between 
about 80^N and 80^S. 

In the routine data processing (Chu, 1986), the radiance data are first 
converted into a slant path transmission. They are then inverted to obtain 
vertical profiles of the aerosol extinction at each wavelength as well as the 
concentration of the gaseous constituents. The first step in the inversion 
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TABLE 1.1. SAGE II. Spectrometer Characteristics 
(McMaster, 1986). 


Primary 

Channel Constituent 
Number Detected 


Center 

Wavelength Bandwidth 

(nm) (nm) 


1 

Aerosol 

1019.7 

19.6 

2 

H 2 O 

935.5 

20.0 

3 


599.9 

14.5 

4 

Aerosol 

525.0 

14.8 

5 

NO 2 , Aerosol 

454.4 

1.9 

6 

CM 

0 

447.5 

3.2 

7 

Neutral Density 5 

384.6 

19.8 


Aerosol 
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procedure is to separate the contribution to the slant path transmission from 
each specie, as well as that due to the molecular atmosphere at each of the 7 
wavelength bands. Following this, the individual vertical profiles for each 
specie is obtained. Separation of contributions from each specie is complex 
and troublesome since for all channels there is a contribution to the optical 
depth from not only the main species but also from the other species. When 
one species has a large signature which extends across a wide wavelength 
range, as for example occurs when a high altitude cloud is present, there may 
be considerable uncertainty associated with the inverted product for other 
species. An estimate of this uncertainty is made along with the inversion 
and forms part of the data product. 

A major part of the total post-launch data validation process is the 
operation of field experiments in which correlative data are collected on the 
aerosols and gaseous constituents. These measurements have been made from 
both ground based and airborne platforms. The former includes lidar 
measurements; the latter includes balloon and aircraft in situ particle 
counters and lidar, rocket and balloon-borne ozone measurements and 
balloon-borne nitrogen dioxide and water vapor measurements. In the present 
study, we have made use of such measurements, as well as data available from 
the SAM II satellite and the historic data set available from the SAGE I 
satellite (McConnick et al. , 1979) . The p resent data validation process also 
involves the screening of the SAGE II data p roduct , both by exami nation of 
tables and plots of the inverted data and by actual data use. A final aspect 
of the data validation is the analysis of associated but not necessarily 
correlative data sets, collected by many different techniques, to provide a 
framework or climatology against which the SAGE II product may be assessed. 
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1,3 Data and Validation Status 


The process of data validation for SAGE II has been to a large extent an 
iterative one. Initially, only a small quantity of SAGE II data was 
available. As a preliminary data product, this contained many errors and 
omissions of varying significance. The more obvious of these were quickly 

found by the authors of this report, and co-workers, resulting in new 

improved versions of the SAGE II data product. The data product on which the 
majority of the work described in this report is based, covers the SAGE II 
measurements from October 1984 (launch) to May 1985 and two correlative 
experiments, one performed in Laramie, Wyoming in November 1984, and a field 
experiment performed in Alaska in August 1985. It should be noted that^th^ 
current dat a set is not final. While, as a result of the feedback process of 
data validation, this data set is vastly improved as compared to the original 
verion, some errors and areas of uncertainty still exist. A new corrected 
version of the data set will be ready in the near future. It is anticipated 
that this will be acceptable for scientific use and that it will be the first 
data set with archival status. 

In this report, we have not attempted to list the many minor errors 
found in the initial data product. Some major errors which still exist, 

e.g*, the sunrise-sunset effect which is still present in the current data 
product, are described. There is also an indepth discussion of analysis of 
some of the correlative data, notably the aerosol size distribution 

measurements and comparison with the SAGE II multiwavelength extinction 
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data. This report also contains the results of scientific analyses of the 
SAGE II data, associated data sets and related theoretical studies* It 
should be emphasized that such studies, based on the current preliminary SAGE 
II data set, are not final* While the qualitative aspects of the results 
obtained are not expected to be in error, quantitative values will be subject 
to change when the analyses are repeated on the revised data set* 

1*4 Tasks and Organization of Report 

The statement of work specifies seven task areas. Work has been carried 
out under each of these areas and this work is reported in sequence in 
Sections 2 through 8 of this report. The title given to each main section is 
the same as the work statement for the corresponding task. Although work has 
been carried out under each task area, the division is not necessarily 
uniform. This division represents the results of discussions between the 
contractors and the NASA Technical Monitor as to the work areas requiring 
major attention. Subsections under each section correspond to separate 
studies within the task area. 

Some of the work carried out has been presented at scientific meetings 
and a small part has been accepted for publication. That the latter is not 
larger is attributed to the fact that to date all current analyses of SAGE II 
data are based, as explained in the previous section, on a preliminary data 
set. As such, it is not yet suitable for formal publication. Where material 
has been presented or published, the corresponding abstract or paper is 
included in the Appendix to the report. To avoid duplication, a brief 
outline is made in the body of the report and reference is made to the 
appropriate Appendix. 
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2 . TASK 1 — mVESTIGATE THE GLOBAL CXIMATOLOGT OF STRATOSPHERIC AEROSOLS, O^, 

WATER VAPOR AND NO^, INCLODING SPATIAL AND TEMPORAL VARIABILITY 

2.1 Stratosj^eric Aerosol Climatology 

As part of the work associated with the validation of SAGE II data, the 
stratospheric global aerosol climatology from 1979-1985 has been reviewed. 
This period has included many volcanic eruptions that have affected the 
stratosphere. The stratospheric aerosol has been intensely studied by 
ground-based, airborne, and satellite techniques with a resulting increase in 
our understanding of its climatology. This work has been summarized and 
presented as an invited review paper at the Symposium on Radiation in the 
Middle Atmosphere at the lAGA-IAMAP meeting in Prague, Czechoslovakia, held 
on August 5-17, 1985. A copy of the abstract for this paper is included in 

this report as Appendix 1 (Kent & McCormick 1985a). 

As part of the same study, the aerosol climatology, as shown by SAGE II, 
has also been examined for compatibility with the results established from 
SAM II, SAGE I, and airborne lidar. Some of the more detailed results of 
this study will be shown elsewhere (Sections 2.2 and 7.1); here, it may be 
noted that the preliminary SAGE II data shows many expected characteristics. 
The most prominent among these are: 

1. the high aerosol concentrations still present in 1984 and 1985, 
following the eruption of El Chichon in 1982, 

2. a decrease in concentration between October 1984 and May 1985 (Limits 
of SAGE II data currently available for study) . 
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3, A zonal structure to the aerosol concentration similar to that seen 


after earlier volcanic eruptions (Kent and McCormick, 1984)* 

2.2 Software for SAGE II Glc^l Plots 

As part of work carried out under a separate contract (NASI- 17032 ) , 
software has been developed to enable global maps of SAGE I 1 pro optical 
depth to be produced* At first sight, this appears to be a relatively simple 
matter as the satellite data systematically sweeps in location from one 
extreme of latitude to another. In practice, difficulties arise because of 
the following; 

1. The satellite data points do not fall on a regular grid. Although 
the longitudinal interval is fairly regular ( -^24^ intervals), the latitude 
interval varies greatly between high and low latitudes. Moreover, there is a 
longitudinal shift in data location on successive days. 

2. Missing data are not infrequent, either requiring an interpolation 
routine or resulting in areas of missing data on the final map. The latter 
almost certainly means that contouring routines cannot be used on the map. 

The SAGE I software has been converted for use with the SAGE II data 
product. The purpose of this conversion was (1) to obtain a useful tool for 
scientific investigation of global aerosol climatology, and (2) to obtain a 
display which permits rapid detection and screening of anomalies in the data 
itself. The algorithm for production of these maps involves the following 
steps . 


2-2 


K Selection of a suitable portion of SAGE II data for mapping, e.g., a 
latitudinal sweep between chosen dates — this should be either a sunrise or a 
sunset sequence. 

2. Production of the variable to be plotted at each data location (in 
the present case, 1 pm optical depth). 

3. Determination of the limits of the latitude sweep and the locations 
of any missing data within that sweep. 

4. Interpolation of missing data, provided the data gaps do not exceed a 
certain predetermined size. The interpolation involves the data points 
neighboring the data gaps. 

5. Interpolation of the filled data sequence to a x global 

latitude-longitude grid. This is a fairly complex procedure involving 
definition of quadrilateral data cells within the SAGE II data product, 
followed by a projection onto a retangular grid and finally a standard 
interpolation to the unit cell structure. 

6. Use of the 1^ x 1^ data matrix to produce either a shaded or a 
contoured data product, or one containing both. The data matrix may also be 
used to produce color-coded maps. 

User friendly control procedures have been written for use with the 
above programming sequence. Although these are, at present, confined to the 
production of 1 pm optical depth maps, their extension to other data products 
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appears to be a relatively simple matter. This could include, for example, 
optical depth and extinction at other wavelengths, or water vapor and ozone 
concentrations. The procedure then becomes a very useful and powerful tool 
for displaying and screening the data product. 

An example of the data product is shown in Figure 2.1. This figure shows 
the I pm optical depth, measured from 2 km above the tropopause, for the 
period October 24 - November 20, 1984. The data shown are for sunset events 

only . 


Six levels of shading are available and contours have been added to 
delineate the shaded regions. The high latitude resolution and quality of 
the data product are apparent. In scientific terms, this figure shows the 
very strong banded structure apparent in the stratospheric aerosol. At this 
time, the aerosol is a well-dispersed and partially decayed aerosol derived 
mainly from the El Chichon eruptions of April/May, 1982. The zonal structure 
has been observed previously both in early post-volcanic aerosol (during 
1981) and in pre-volcanic aerosol (during 1979). 

2.3 Analyses of Error Bars on 1 pm Aerosol Data 

An important SAGE II data product is the error bar which is supplied for 
each data point. Much subsequent scientific analysis will involve screening 
the data using not only the quality factor which is represent at ive of a 
complete vertical profile, but also the individual error bars for each 
altitude . As part of our examinat ion of the preliminary data set , an 
analysis was made of the error bars associated with the SAGE II 1 pm 
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Figure 2.1. SAGE II measurements of 1 jim optical depth from 2 km above the 
tropopause, October 24 - November 20, 1984. 
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extinction data. Two aspects of the data were examined: (1) Did the error 
bars appear to give a reasonable representation of the statistical 
fluctuations in the data set and (2) How did the SAGE II error bars compare 
with those given in the SAGE I and SAM II data sets ? Programs were written 
for all three satellite data sets which (I) divided the measured data in 10^ 
latitude bands and (2) calculated the mean 1 pm extinction and extinction 
ratios within each band, the root mean square variation of these quantities 
within the band, and the mean error bar. Assuming little or no systematic 
error exists within the data set, the error bar is expected to be less than, 
or approximately equal to, the rms variation within the data set. This 
follows from the fact that the latter includes both data errors and actual 
geophysical variations. The expected inequality is not found to be always 
true for either the SAGE I or the SAGE II data set. For the SAM II data set, 
which is inverted by a different scheme, the rule is generally true. 

Figure 2.2 shows typical results for SAGE I and SAGE II taken for 
identical months (in different years) and for the same latitude band. The 
extinction ratio profiles for the two periods are significantly different 
below an altitude of 25 km, due to the presence of aerosol from the eruption 
of El Chichon in 1982, which produced an increased extinction ratio as shown 
in Figure 2.2(b). There is also a significant difference in the extinction 
ratios above 25 km which is not due to the volcano. Part of this difference 
may be due to the use of a slightly incorrect molecular atmosphere for SAGE 
I. (W. Chu, 1986 private communication.) 
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Figure 2.2. Altitude variation of mean 1 |im extinction ratio, extinction 
ratio error and RMS variation within a 10° latitude band for 
(a) SAGE I data, April 1979; (b) SAGE 11 data, April 1985. 
Both data sets are for latitude bands 40° — 50°N. 
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Below 25 km, the mean error bar for both SAGE I and SAGE II data is of 
the same order or less than the observed rms variation within the data set 
and no serious discrepancy exists. Above 25 km, the error bars increase 

rapidly for both satellite data sets; in contrast, the rms variation 
decreases. A major part of this behavior difference is due to the 
introduction of vertical smoothing into the data inversion process when the 
extinction falls below a certain predetermined value. This smoothing is 
carried out over a vertical height interval of 5 km; the error bars are not 
smoothed and would thus be expected to be larger (relative to those at a 

lower altitude) by a factor of sfs* This factor is probably sufficient to 
account for the relative growth in the error bars (for both SAGE I and SAGE 
II) up to an altitude of about 30 km. Above this altitude the error bars 
appear to be an overgenerous estimate of the errors within the data set. In 

the case of the SAM II where a different inversion scheme was used, no 

comparable height variation in the significance of the error bar was found. 
This vertical variation in the significance of the SAGE I and SAGE II error 
bars indicates that care must be taken in any analysis that depends upon 

them, particularly as it relates to the altitude range 20-30 km. 

2.4 Ozone and NO 2 Diurnal Variations 

In this section, we report the results of validating the available 

preliminary SAGE II data set in terms of examining the diurnal variability of 
the mesospheric 0^ and stratospheric N02* There are three occasions in the 
first eight-month sampling period of the SAGE II instrument at which the 

sunrise and sunset measurements cross each other at certain latitudes. 
Figure 2.3 shows the location and the time of the occurrence of these 
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occasions (the circled numbers). Figure 2.4 shows the averaged profile for 
the sunrise events (solid line), and also that of the sunset events (dashed 
line) for such an occasion which happened in May 1985. Only those sunrise and 
sunset measurements for which the differences in latitude are less than 5^ 
are included in the analysis. As indicated in Figure 2.4, the averaged 
sunrise profiles show about a 6% higher 0^ concentration than that of the 
sunset profiles in the altitude range between the lower limit of the figure 
(30 km) and approximately 60 km. Since in the stratosphere the 0^ diurnal 
variation is negligibly small (e.g., Wang et al., 1980), Figure 2.4 indicates 
a bias in the SAGE II 0^ data set. 

The corresponding averaged sunrise and sunset profiles for NO^ are 
displayed in Figure 2.5. The cross behavior of the NO^ profiles at altitudes 
-25 km and 37 km is less anticipated and is difficult to interpret. In 
general, the stratospheric photochemistry suggests a similar vertical 
distribution of NO 2 concentration with the sunset NO^ higher than that of 
sunrise (Chu and McCormick, 1986). As a result, further investigation on the 
preliminary SAGE II NO^ data set is required. 

Because of the unexpected behavior difference between the sunrise and 
sunset 0^ and NO 2 , derived from preliminary SAGE II data set, this type of 
data comparison has been extended to all the other SAGE II measurements, 
including H^O and the four aerosol channels. The results of this comparison 
are given in Figures 2.6a to e, respectively. As can be seen, the results 
also indicate a systematic bias with the sunrise measurements at a higher 
altitude than the corresponding sunset values. 
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Figure 2. 4, Comparison of preliminary SAGE II ozone prof iles for sunrise and sunset measurements. 
Data are for May 1985 and averaged over several events. 



o 



m m in 

m <t ^ 


in m ^ 


(W)i)aani.inv 


CO (U 


2-12 


Figure 2.5. Comparison of preliminary SAGE II NO 2 profiles for sunrise and 
measurements. Data are for May 1985 and averaged over several 
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Figure 2.6(a). Comparison of preliminary SAGE II H 2 O profiles for sunrise and sunset 
measurements. Data are for May 1985 and averaged over several events 
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Figure 2.6(b). Comparison of preliminary SAGE II 1.02 pm aerosol extinction profiles 
sunrise and sunset measurements. Data are for May 1985 and averaged i 
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Figure 2.6(c). Comparison of preliminary SAGE II 0.525 pm aerosol extinction profiles 

for sunrise and sunset measurements. Data are for May 1985 and averaged 







In order to further confirm the above findings in this preliminary May, 
1985 SAGE II data set, the profiles obtained at the other two adjacent 
sunrise and sunset crossings, as shown in Figure 2.3, have also been 
examined. Again, a similar unexpected behavior of the sunrise and sunset 
profiles of this preliminary SAGE II measurement is shown in the averaged 
profiles for these two occasions* 

Relevant to the above findings of the bias in the current SAGE II data 
set, G* Yue of NASA-LaRC and M* T. Osborn of SASC Technologies Inc. have 

found, independently, that good agreement exists between the 1 pm aerosol 
extinction profiles of individual SAGE II sunset events and SAM II events 
that occurred close together in space and time. This is not true for 

similarly collected SAGE II sunrise and SAM II observations. In their 
analyses , the SAGE II sunrise 1 pm aerosol extinction profiles are 

systematically shifted to a higher altitude as compared to those for SAM II, 
Thus , from their results , and that of the sunrise and sunset analysis 
mentioned earlier, it is most likely that the bias exists in the preliminary 
SAGE II sunrise data set. 

2,5 Anoaalous Values of the 0,525/1 #0 pm Extinction ratio 

As part of the development of an aerosol climatology for the 

stratosphere and upper troposphere , a study has been made of the wavelength 
extinction ratio (Extinction (X= 0,525 pm) /Extinction ( X= 1,02 pm)). This 
ratio is related to the particle size distribution (see Section 3,1.1,). For 
small particles such as exist in the statosphere, a ratio in the range 2-5 is 
to be expected. For cloud particles, which are much larger, the ratio should 
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be about unity. It is very difficult to construct a realistic particle size 
distribution for which this ratio is substantially less than unity, although 
values slightly less than unity are possible (Yue, Private Communication, 
1986). In the current study, SAGE II aerosol data have been divided into 20^ 
latitude bands and then binned, according to altitude from 0 to 30 km, and 
according to 1 pm extinction ratio from 0.01 and 1000.0. The data have been 
filtered to remove unreliable data points using both the quality factor 
supplied with the inverted profiles, as well as individual error bars at each 
altitude. Most of the 0.525 pro/1.02 pm extinction ratio values are as 
expected, but some anomalies have appeared. The results may be summarized as 
follows; 

1. In the stratosphere, most wavelength extinction ratio values lie 
between 2 and 6; the higher values, associated with smaller particles, occur 
at the greatest altitudes. 

2. In the free troposphere, more than one aerosol mode is apparent. 

a. Cloud Mode » Wavelength extinction ratio values are close to unity 
as expected. 

b. Coarse Aerosol Mode (Northern hemisphere) . Wavelength extinction 
ratio values are about 2. 

c. Fine Aerosol Mode (Southern hemisphere) . Wavelength extinction 
ratio values are in the range 3-4. 

In addition, however, many values substantially less than unity are 
found (as low as 0.1). These cannot be true values (private communication 
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with W, P. Chu). These values probably arise during the inversion process. 
They only occur at the lower altitudes and are associated with the occurrence 
of high extinction or optical depth values* Their origin is being examined 
and it is anticipated that they will be eliminated from the next version of 
the SAGE II data set* 
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3. TASK 2—INVESTIGATE THE CORRELATION OF STRATOSPHERIC AEROSOLS, WATER 
VAPOR AND NO^ WITH EACH OTHER AND VARIOUS METEOROLOGICAL PARAMETERS 

3.1 Modeling the Stratospheric Aerosol Extinctions at 1*0 and 0.45 pm 

3. 1. 1 Introduction 


Since satellite instruments, such as SAGE II, provide aerosol extinction 
information at several wavelengths, model calculations are useful in 
interpreting the observed results. Furthermore, model analyses are also 

essential to the retrieval of aerosol properties from these multiwavelength 
extinction measurements (e.g., Yue and Deepak, 1983; see also Section 6.1). 
More specifically, one needs to know clearly the dependence of the 
multiwavelength extinction on the aerosol properties including aerosol size 
distribution, number concentration, refractive index, etc. 


In the stratosphere, several different models of the aerosol size 
distribution have been introduced (e.g., Russell et al., 1981). Perhaps the 
most widely used one is the log-normal representation introduced by Pinnick 
et al. (1976) which is given by ^ ^ 


n(r) = 


/27T £na 


exp 


£n^(r/rg) 

2£n^a 


(3.1) 


where n(r) is the number of particles per cm for particles whose radii are 

-3 

between r and r + dr, N is the total number concentrations (cm ), o is the 

o 

width of the log normal curve, and r is the mode radius. The values 1.86 
^ ’ 8 

and 0.0725 pm have been used commonly in the literature for a and r , 

6 
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respectively. The aerosol refractive index, on the other hand, depends on 
the composition of the sulfate droplets (Gmitro and Vermeulen, 1964). 
Furthermore, it is also temperature dependent, following the Lorentz-Lorenz 
relationship (Longhurst, 1964). For example, the value 1.43 corresponds to 
the composition of 75% of H^SO^ and 25% of H^O at about 218^K (Pinnick et 
al., 1976). 

For stratospheric sulfate droplets, the particles will absorb (or 
evaporate) water rapidly in response to changes in the ambient temperature 
such that the equilibrium will be reestablished between the environmental 
water vapor and the vapor pressure under the new composition of the sulfate 
droplets. The net result is a change in both the composition and the size of 
the particles (e.g., Steele and Hamill, 1981). 


Recently, Steele and Hamill (1981) have developed a theoretical model 
for determining the effect of changes in temperature and water vapor on the 
aerosol particle size. In addition, this model has been applied in 
conjunction with the Mie scattering theory to determine the dependence of the 
aerosol extinction at 1.0 pm and the aerosol backscatter at 0.6943 pm on the 
temperature and humidity in the stratosphere. The log-normal size 


distribution given by Eq. (3.1) is used in the calculation of aerosol optical 

-3 

properties with = 10 cm , a = 1.86 and r^ = 0.0725 pm. They have further 
extended their model to investigate the formation of Polar Stratospheric 


Clouds (PSCs) based on SAM II 1.0 pm extinction observations (Steele et al.. 


1984). 


On the other hand, Yue and Deepak (1983) have developed an effective 
method for retrieval of stratospheric aerosol size from satellite extinction 
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measurements at two wavelengths, e.g., 1.0 pm and 0.45 pm. This method has 

been applied to the aerosol extinction data set derived from the SAGE I 
measurements (Yue and Deepak, 1984). The development of their method is 

mainly based on some observed stratospheric features. More specifically, it 
is assumed that the overall stratospheric temperature and the aerosol 
particle composition (Wg, weight percentage of H 2 S 0 ^) and the width of the 
log-normal curve are 200^K, 75%, and 1.86, respectively. It should be noted 
that , with the above specified values of these relevant parameters , the 

interpretation of the different retrieved aerosol sizes would be due to the 
different gaseous concentrations of ambient water vapor and the mass content 
of sulphuric acid in the droplet according to the theoretical model of Steele 
and Hamill (1981). It is important to keep in mind that, in their model, the 
aerosol droplets are assumed to be in equilibrium with respect to water 
vapor. The condition for droplets to be in equilibrium with respect to 

sulphuric acid is not necessary. This is because the t ime taken for the 

droplet to undergo a mass increase of 1/e by water absorption is of the order 
of seconds, whereas it is at least 10^ times larger for the sulphuric acid 
absorption. Thus, the growth by sulphuric acid absorption becomes important 
only for processes occurring over long time periods. Furthermore, since the 
vapor pressure of water over a flat solut ion of a g iven solute weight 
percentage (W^,) is a function of temperature only, and since the Kelvin 
effect in Che vapor pressure is important only for very small aerosol 
particles, all the droplets would be characterized by a specific if they 
are at the same ambient temperature and water vapor concentration. This 
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implies that a large particle would contain more sulphuric acid for a given 
ambient water vapor concentration (W). 

As mentioned earlier, Steele and Hamill (1981) have used a= 1.86 and r 

g 

= 0.0725 jjm in their calculation of the aeorsol extinction at 1 pm. These 

constants are also introduced into their study of the formation of Polar 

Stratospheric Clouds. Yue and Deepak (1984) have assumed T = 200 K, W = 75%, 

and a = 1.86 in their development of the aerosol size retrieval method from 

measurements at two wavelengths. In reality, however, the stratospheric 

temperature may have a range from 185^K to 230^K. Similarly, the values of a 

and r may change from place to place and time to time. For example, a mode 
g 

radius r of 0.097 pm and a spread a of 2.02 are considered to be appropriate 
g 

for the stratosphere following a moderate volcanic event (Chu and McCormick, 
1979). Furthermore, in the free troposphere, a value for O of 2.985 is 
considered appropriate for background aerosol (McClatchey et al., 1980). 

In this section, we shall examine the variations of aerosol extinctions 
at I .0 pm and 0.45 pm for different values of a and at various ambient 
temperatures and water vapor. In addition, this study also examines the 
dependence of the aerosol size on the extinction ratio (0.45 pm to I.O pm) 
under various conditions of atmospheric temperature, water vapor, and the 
width of the aerosol size distribution (a) values. It serves as a 
sensitivity analysis for retrieval studies of the stratospheric aerosol size 
distributions using SAGE II multiwavelength aerosol extinctions. As shown in 
Section 6.1, the retrieved aerosol size distributions from SAGE II 
measurements can be further used for SAGE II data validation purposes. 
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3. 1.2 Coniputation and results 


Since the intention of this analysis is to examine the sensitivities of 
the calculated aerosol extinction at I.O pm and 0.45 pm wavelengths, at given 
temperatures, with respect to changes in water vapor concentration and 
aerosol size distribution, consideration of the influence of dynamics is not 
necessary. In the analysis, the theoretical model developed by Steele and 
Hamill (1981) is employed for computation. The log-normal representation of 
Eq. 3.1 is used to determine the aerosol optical properties according to the 
Mie scattering theory. In addition, we discuss only the case of supercooled 
droplets at sufficient low temperature, although Steele and Hamill (1981) 
have concluded that, in the case of PSCs, the particles are more likely to be 
frozen ones. 

Figures 3.1a to 3.1c show the calculated results of aerosol extinction 

at 1.0 pm for fixed 0 values of 1.674 , 1.86, and 2.064, respectively. The 

value 1.86 is the most widely used value for describing the background 

stratospheric aerosol. The lower and higher values of a differ from this 

value by -10% and +10%, respectively. For each fixed a, the extinction is 

calculated for mode radius (r ) with values 0.05, 0.07, 0.09, 0.11, and 0.13 

8 

pm and water vapor mixing ratios of 2 , 4 , and 6 ppmv . These water vapor 

-4 -4 

concen trat ions correspond to the vapor pressure of2xl0 ,4x10 , and 6 

-4 

X 10 mb , respectively, at the 100 mb atmospheric pressure level . Figures 
3.1a to 3.1c show many interesting features. First of all, the extinction 
for temperatures greater than about 210°K is not sensitive to changes in 
ambient water vapor pressure but is to the mode radius (r ). This condition 
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ppmv at 100 mb for mode radius rg = 0.05, 
0.09, 0.11, and 0.13 pm, with fixed values of 
the lognormal aerosol size distribution: 

- 1.674. 




3 -: 


temperature for water vapor mixing ratios of 2, 4, 
and 6 ppmv at 100 mb for mode radius rg = 0.05, 
0.07, 0.09, 0.11, and 0.13 pm, with fixed values of 
0 in the lognormal aerosol size distribution; 

(b) a = 1.86. 
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d 6 ppmv at 100 mb for mode radius rg = 0.05, 

07, 0.09, 0.11, and 0.13 pm, with fixed values of 
in the lognormal aerosol size distribution: 

) a = 2.064. 



becomes reversed at sufficiently low temperatures as the aerosol particles have 

grown by absorbing large amounts of water molecules so that they become 

almost pure water droplets. These features have been noted by Steele et al, 

(1984). Figure 3.1 also depicts the sensitivity of the extinction to changes 

in a. This sensitivity is better described in Fig. 3.2. Figure 3.2 shows the 

variation of the 1.0 pm extinction with temperature for a fixed mode radius 

r (= 0.0725 pm) at the three constant a values mentioned earlier, i.e., 

8 

1.674, 1.86, and 2.046. For each of the O values, computations are carried 

out again for water vapor concentrations of 2, 4, 6 ppmv. Figure 3.2 shows 
that a 10% change in Q may lead approximately to a factor of 3 change in the 
extinction for temperatures above a certain threshold value, below which the 
droplet becomes a supercooled particle. 

The results of the corresponding analysis for wavelength 0.45 pm are 

given in Figures 3.3(a-c) and 3.4 in terms of the wavelength extinction ratio 
(defined as the ratio of extinction at 0.45 pm to that at I.O pm). One may 
notice from Figure 3.3, that for a fixed a, a doubling of the value of r 

8 

leads to a reduction of this ratio by a factor of approximately 3, where a = 

1.674, for a given water vapor concentration and a given temperature above 

210^K. Figure 3.3 also shows that this reduction decreases as a increases. 

Figure 3.4 is the same as Figure 3.2, except that it is for the wavelength 

extinction ratio. It indicates that the sensitivity of the extinction ratio 

to changes in 0 for fixed values of r , temperature, and water vapor 

8 

concentration depends on the value of 0 itself. A decrease of a from 1.86 to 
1.674 leads to an increase of the ratio by a factor of approximately 2 in the 
case of r = 0.0725 pm. 
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and 6 ppmv at 100 mb for rg = 0.725 pm and O = 1.674 
1,86, and 2.064 in the lognormal aerosol size 
distribution. 






0.45 ]im to that at 1.0 ymi) as a function of temperature 
for water vapor mixing ratios of 2, 4, and 6 ppmv at 
100 mb for mode radius rg = 0.05, 0.07, 0.09, 0.11, and 
0.13 pm with fixed values of a in the lognormal aerosol 
size distribution: (b) a = 1.86. 
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for water vapor mixing ratios of 2, 4, and 6 ppmv at 
100 mb for mode radius rg = 0.05, 0.07, 0.09, 0.11, 
0.13 pm with fixed values of o in the lognormal aero 
size distribution: (c) a = 2.064. 



250 



Figure 3.4. The aerosol extinction ratio (aerosol extinction at 

0.45 pm to that at 1.0 pm) as a function of temperature 
for water vapor mixing ratios of 2, 4, and 6 ppmv at 
100 mb for rg = 0.0725 pm and a = 1.674, 1.86, and 
2.064 in the lognormal aerosol size distribution. 



To further examine the relationship between the mode radius r , and the 

8 

wavelength extinction ratio under various values of water vapor 
concentration, o , and temperature, the aerosol extinction ratio has been 
determined for a range of values of r from 0.01 pm to 0.30 pm. The results 

s 

of this calculation for a temperature 190^K are shown in Figure 3.5. Like the 

study given in Figure 3.2, the analysis involves three values of CT , i.e., 

1.674, 1.86, and 2.046. For each of these values, the sensitivity study is 

carried out with respect to four different values of water vapor 

concentrations (from 2 ppmv to 8 ppmv) . In Fig. 3.5, the results for 2 ppmv 

are given by the solid lines. The dashed lines correspond to 4 ppmv, and the 

solid-dashed lines are for 6 ppmv and 8 ppmv. First of all, the extinction 

ratio shows a monotonic increase with the decrease of the mode radius. 

Similar features can be found in the analysis by Yue and Deepak (1983, Figure 

3.1a). Figure 3.5 also shows that the extinction ratio is sensitive to the 

value of a, especially when the mode radius is small (e.g., less than 0.1 

pm). As to the influence of the water vapor concentration, the ratio does 

not seem to be sensitive to its variation in the range from 2 ppmv to 8 

ppmv. The reason that the two higher water vapor concentration cases, in 

particular, show very similar values for the extinction ratios is that, for 

the given r , temperature (190^K) and a, the aerosol particles are 
8 

supercooled droplets in this model calculation. As a result, they have very 
similar optical properties. From Figure 3.5, it is clear that the extinction 
ratio is rather insensitive to the water vapor concentration in the range 
from 2 ppmv to 8 ppmv at a given temperature. Similar computations have been 
made for temperature at 200^K and 210^K. The results of these computations 
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Figure 3.5. The aerosol extinction ratio (aerosol extinction at 
0.45 um to that at 1.0 pm) as a function of mode 
radius rg for G = 1.674, 1.86, and 2.064 in the log- 
normal aerosol size distribution. The temeprature i 
fixed at 190°K. 



differ only slightly from those shown in Fig, 3,5, In other words, the 
extinction ratio is not very sensitive to changes in temperature, 

3. 1.3 Discussion. 


Satellite technology gives us an effective means of studying 

stratospheric aerosols with global scale coverage. These measurements are 

accomplished through de terminat ion of the opt ical properties of aerosol 

particles . The interpretation of the measured optical properties , as well as 

the retrieval from satellite extinction measurements of the essential 

parameters of these particles, such as size distribution, composition etc., 

rely on detailed model analyses . In this section, we have examined the 

effect of temperature, water vapor, and aerosol size distribution on the 

aerosol extinctions at 0.45 pm and 1.0 pm wavelengths, and their ratio (0.45 

pm to 1.0 pm). It is concluded that the extinctions at these two wavelengths 

are not sensitive to changes in ambient water vapor pressure but to the mode 

radius (r ) for temperatures greater than about 2I0^K. This condition becomes 

reversed at sufficiently low temperatures. If we vary the width of the size 

distribution a , it is found that a 10% change in a may lead approximately to 

a factor of 3 change in the extinctions. As to dependence of the wavelength 

extinction ratio on ambient temperature , water vapor concent rat ion , and the 

size distribution, it is shown that, for sufficiently small mode radius (r < 

g 

0.1 pm), the extinction ratio is sensitive to changes in mode radius. This 
sensitivity increases with the decrease of the mode radius. It is also found 


the wavelength ext met ion ratio is 


amuieiiu 


temperature and water vapor concent ration but to the width of the size 
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distribution (a), especially when the mode radius is less than O.I pm. As a 
result of this sensitivity, one can only derive an equivalent mode radius 
from the extinction information at two wavelengths (i.e., 0.45 pm and 1.0 
pm). This feature has been noted in the study by Lenoble et al. (1984). 
They mentioned that, roughly, an increase of O from 1,6 to 1.92 gives results 
equivalent to an approximate reduction of radius (r ) by a factor of 2. 
Obviously, to resolve this problem, it is necessary to incorporate at least a 
third independent extinction measurement. With this additional information, 
it is possible to determine a appeared in Eq . 3.1. As shown in Section 6.1, 
the SAGE II multiwavelength aerosol extinction measurements allow us to 
derive q well as r in the retrieval analysis. Furthermore, the 

s 

measurement from the SAGE II water vapor channel will also lead to 

information about the composition of the stratospheric aerosols. 

3.2 Dependence of Aerosol Composition on Water Vapor and Temperature 

The preliminary results of the SAGE II water vapor measurements during 
April 1985 are being used in conjunction with the temperature profiles 
provided by NOAA to investigate the chemical composition of stratospheric 

aerosols . This investigation also serves as an additional data validation 
analysis for the SAGE II water vapor measurement . In this analysis , the 

water vapor partial pressures at a given standard pressure level in the lower 
stratosphere are obtained from the SAGE II measurements using an 

interpolation scheme. The interpolated results are set in a 2-dimensional 
space , (i.e., water vapor pressure against temperature) . The results of a 

typical example of the analysis at 100 mb pressure level are given in Figure 
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3.6. Figure 3.6 shows also the distribution curves of equilibrium water vapor 
pressure as a function of temperature for pure water, pure ice, and six 
different H^SO^ weight percentages of sulphuric acid-water solutions. As one 
can see, the majority of the April observations show an equilibrium water 
vapor pressure of about 3 x 10 mb, and indicate a weight percentage of 

70Z or so for aerosol droplets. In Figure 3.6, several SAGE II measurements 
of water vapor pressure show values lower than 3 x 10 ^mb. These values are 
unlikely to be real for the reason that the associated air parcels would have 
to go through an environment with temperatures below ISO^K to have most of 
the water removed. These low SAGE II water vapor results may have occurred 
due to over-correction for the aerosol effects on the water vapor channel 
measurements (private communication with W. P. Chu) . 

Figure 3.6 also indicates that many SAGE II water vapor measurements 
suggest relatively low values of percentage of the aerosol 

droplets. These aerosol particles may be of tropospheric origin and are 

located in the tropics. In order to see the vertical variations of the 
aerosol composition, the preliminary SAGE II water vapor information and the 
associated temperature profile have been used to estimate the percentage 
weight of H2S0^ in the aerosol droplets. We have plotted each vertical 
profile sampled in the tropical region separately from that of northern high 
latitudes (> 50^N) to gain some understanding of the latitudinal variation. 
These results are displayed in Figures 3.7a and 3.7b, respectively. As shown 
in Figure 3.7a, all the typical profiles are featured generally by a 

down-slope from the upper right of the figure, corresponding to high 
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Figure 3.6. The vapor pressure as a function of temperature over water, ice, and six 
different H 2 SO 4 weight percentages of sulphuric acid-water solutions. 

The crosses correspond to the SAGE II water vapor measurements and the 
associated temperatures. 
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Figure 3.7a. Distributions of the preliminary SAGE II water vapor measurements between ± 10® 
latitude and associated temperature (provided from NOAA) , April 1985 in the 
2-dimensional space of temperature and vapor pressure. Shown also in the figure 
are the saturation vapor pressure of pure water, ice, and six different percentages 
of H^SO, (by weight) of sulphuric acid-water solutions. 




Figure 3.7b. Same 











temperature and high pressure, to the low temperatures 200^K) as the 

altitude increases. As the altitude increases futher, all the profiles 
become merged together at a location specified approximately by T«2I0^K and 

-4 

V»2.5xl0 mb. These profiles then slope down toward the right as the altitude 

continues to increase. It must be pointed out that the part of the profile 

which shows higher water vapor pressure than that for the pure water is 

unlikely to be real. Similarly, the portion of the data with the water vapor 
— 4 o 

pressure below ~10 mb near T-200 K is probably questionable. 

In the case of April 1985 SAGE II water vapor at northern high latitudes 
(> 50°N), Figure 3.7b indicates generally that the temperature and water 
vapor pressure decreases first as the altitude increases until a temperature 
*'200°K is reached. This characteristic temperature must correspond to the 
tropopause temperature at this high latitude region. Since the altitude of 
this temperature extends up several kilometers, the vertical profiles of the 
water vapor on T-V diagram show a drop of water vapor pressure with 
relatively less temperature changes as the altitude increases in this 
particular altitude region. The further increase in altitude results in 
further decreases in water vapor pressure and increases in temperature. 
Thus, the profiles on the T-V plot show a sloping down towards the lower 
right corner similar to that in the tropics (Figure 3.7a). 

In conclusion, we have found that the behavior of the water vapor and 
temperature in the tropics for the month of April is generally different from 
that at high latitudes. In addition, the existing SAGE II preliminary water 
vapor data set at high latitudes looks more reasonable than that in the 
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tropics. However, since the lower water vapor values may be the result of an 
aerosol correction problem and since the aerosol loading from El Chichon 
seems to be decaying (see Sect ion 4.1), the SAGE II water vapor data should 
become more reasonable with time than exhibited in the April 1985 data set. 
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4. TASK 3 — INVESTIGATE THE INJECTION OF VOLCANIC MATERIAL INTO THE 

STRATOSPHERE INCLUDING GLOBAL LOADING AND TRANSPORT 

4.1 Stratospheric Effects of the El Chichon Volcanic Eruption 

The most significant event affecting the stratospheric aerosol that has 
occurred within the past twenty years was the series of eruptions of the El 
Chichon volcano in March/April 1982. These eruptions unfortunately occurred 
after the failure of the SAGE I Application Explorer Mission Satellite power 
supply and before the launch of SAGE II. The stratospheric injection by these 
eruptions has been such as to dominate the stratospheric aerosol for several 
years after the eruptions. It is only recently that the aerosol level, still 
several times that existing at its lowest value in 1979, has dropped enough 
for the effects of other smaller eruptions to be detectable (e.g., that of 
Ruiz in November, 1985). In spite of the fact that global satellite 

extinction data were not available for the earlier phases of the El Chichon 

period, the effects have been well studied. Methods of study have included 
use of data gathered by the SAM II satellite and NASA-LaRC airborne lidar 
system, as well as by numerous other experimental techniques. In view of the 
importance of these eruptions and the relatively strong effects still 
remaining at the time of launch of SAGE II, a review has been prepared on the 
observations of the effects between April 1982 to July 1985. This review has 
included data from SAM II, the NASA-LaRC airborne lidar, as well as early 

data from SAGE II. Where appropriate and necessary, data from other sensor 

systems have been integrated. The review was presented as an invited paper 
at the Workshop on Ozone Variations on Climatological Time Scales: 
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Observations and Theory and Perturbations by UV Flux Variations and the 
Eruption of El Chichon, held in Salzburg, Austria, August 19-22, 1985. A copy 
of the abstract of the review paper is attached to this report as Appendix 2. 
The paper is scheduled to be published as part of a Middle Atmosphere Program 
(MAP) Handbook (Kent and McCormick, 1985b). 
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5« "CASK 4 — XNCLDDIHG THE STUDY OF PHYSICMi MB> CHEMICAL PROCESSES 

5«1 Analysis o£ SAGE 11 and Correlative Measurenents on August 7, 1985 at 

Fairbanks, Alaska 

The presence of aerosols in the lower stratosphere and the upper 

troposphere affects the retrieval of SAGE II water vapor measurements* This 

is also particularly true when there are occurrences of subvisible or visible 

cirrus clouds (Woodbury and McCormick, 1983, 1986)* It is found that during 

the SAGE II ground correlative measurements conducted on August 7, 1985 at 

Fairbanks, Alaska, the aerosol extinction measurements from the SAGE II 

instrument revealed the presence of clouds. This situation is illustrated in 

Figure 5*1. In this section, we discuss the analysis of the formation and 

evolution of aerosols including their physical and chemical properties using 

this SAGE II and correlative data set* Since the satellite instrument is 

viewing the atmosphere from space, it is possible to locate the cloud top 

height within - 0*5km which is the accuracy of the vertical resolution of the 

SAGE II instrument. Figure 5*1 indicates clearly the cloud top is located 

somewhere between altitudes 11 and 12 km* This result is further confimred 

by the correlative dustsonde observations. Figure 5.2 shows the vertical 

profiles of the dustsonde measurements for the N and N channels* The 

1*2 1*8 

distinct peaks of these two channels, which appeared at altitude 12 km, are 
remarkable. These two profiles further indicate that the vertical extent of 
the subvisible cloud is a little over 1 km* 
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Figure 5.1. Preliminary SAGE II aerosol extinction profiles 

obtained on August 7, 1985, at Fairbanks, Alaska. 
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The associated temperature profile for the dustsonde measurements is 
displayed in Figure 5*3 by the dots. Figure 5.3 also shows the temperature 
profile obtained by Oltmans of NOAA/ARL, and that of the SAGE II preliminary 
data set provided by NOAA. Although there is some difference in temperature 
between the one associated with dustsonde measurements and that obtained by 
Oltmans, they both indicate a layer of coldest air located approximately at 
an altitude of 12 km. They also suggest that the vertical extent of this 
coldest layer is about 1 km. Since this is at summer season, it is possible 
that the convective activity leads to upward air motion, and the formation of 
clouds through adiabatic cooling processes. It is also possible that this 
cloud is developed by associating with a frontal system. In order to see the 
details of the development, we have examined the available correlative water 
vapor profiles. The water vapor pressure and the corresponding temperature 
from Oltmans* measurements and that provided on the SAGE II profile tape is 
plotted on a T-V diagram (i.e., temperature against vapor pressure). This is 
given in Figure 5.4. The advantage of using a T-V diagram is that it provides 
an estimate of the weight percentage of stratospheric aerosol 

droplets, and also indicates possible formation of water or ice clouds. In 
Figure 5.4, the five numbers in the parentheses by the symbols indicate the 
altitudes in kilometers of Oltmans* measurements. The location of the 

corresponding data point of the SAGE II observations at the same altitude can 
be identified by the corresponding symbol in this T-V diagram. 

In Figure 5.4, the results of Oltmans* in situ measurements are 
different in two aspects from those Inferred from the SAGE II data set. In 
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VRPOR PRESSURE, MB 



TEMPERATURE, DEG 

Figure 5.4. The vapor pressure as a function of temperature over water. 

ice and six different H 2 SO 4 weight percentages of sulphuric- 
water solutions. The two profiles with symbols are obtained 
from Oltmans’ observations and from SAGE II water vapor and 
NOAA temperature data. The numbers in the brackets indicate 
the altitude in km of associated data points (symbols) . 






the first place, NOAA satellite derived temperature data on SAGE II profile 
tapes is about 5^C higher than Oltmans’ balloon data. Secondly, the 
preliminary SAGE II water vapor data below altitude 15 km are about a factor 
2 greater than Oltmans*, For this reason, the following discussion is based 
on Oilmans* results. It is interesting to note that his data at 10 km are 
located right on the saturation vapor pressure curve for pure ice. This 
feature indicates that the bottom of the cloud layer is about at this level. 
At altitude 11 km, the data point is still very close to the pure ice curve. 
As the altitude increases from 11 km, the observed vapor pressure departs 
from the pure ice curve. The data at 12 km, as one can see, correspond 
approximately to the level of the coldest temperature (also shown in Figure 
5.3). Figure 5.4 further indicates that at 12 km, the aerosol droplets are 
possibly composed of 40% of and 60% of by weight. At 13 km, on the 

other hand, the corresponding composition is about 70% h^SO^ and 30% of h^O. 
In view of the above discussion, based on Figure 5.4, the cloud layer appears 
to extend approximately from 10 km to 11 km. This cloud top location is 
consistent with the results of SAGE II aerosol extinction data set within the 
accuracy of the SAGE II vertical resolution. 

In conclusion, both the SAGE II satellite instrument and the correlative 
measurements, including that made by the dustsonde and Oltmans* water vapor 
sounding, observed the cloud layer during this ground truth event. This 
cloud layer appears to consist of ice particles generated by adiabatic 
cooling processes associated with the upward air motion. Further detailed 
analyses on this cloud development and studies of the aerosol corrections to 
the SAGE II water vapor measurements are recommended. 
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6. TASK 5 — INVESTIGATE THE DETERMINATION OF STRATOSPHERIC MOLECULAR AND 


AEROSOL OPTICAL PROPERTIES USING MULTIWAVELENGTH EXTINCTION AND 

BACKSCATTER DATA AND TECHNIQUES 

6.1 Aerosol Data Validation in Tens of Retrieved Aerosol Size Distribution 
fron SAGE II Multi wavelength Aerosol Extinctions 

As mentioned in Section 3.1, model aerosol size distributions can be 

determined from SAGE II multi wave length aerosol extinctions. To determine 

the model parameters, we have used a nonlinear least squares approach. 
Single mode lognormal as well as modified-ganuna functions have been used in 
the size distribution analysis. In addition, a bi-lognormal size 

distribution can be determined by using both multiwavelength SAGE II aerosol 
extinction data and lidar backscatter data at ruby and green wavelengths • 
The retrieved size distributions can then be used to determine the optical 
properties of the aerosol particles and to calculate the total number density 
of aerosol droplets corresponding to dustsonde measurements. These derived 
results from retrieved size distributions can be used for data validation by 
comparison with correlative measurements. In the following analysis, we 

discuss first the information content of SAGE II aerosol extinctions with 
respect to aerosol size distribution. This is important, since it helps us 
to understand the quality of the retrieved aerosol size distribution. With 
this information, we then discuss the data validation for SAGE II 
observations on two particular days, when correlative measurements were 
available. Part of this analysis in this section has been presented in the 
sixth conference of Atmospheric Radiation at Williamsburg, Virginia on May 
13-16, 1986 (Wang et al., 1986)- 
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6. 1. I SAGE II aerosol size distribution information content 


The aerosol size distribution information content associated with a 

given optical measurement can be described approximately in terms of the 

corresponding efficiency factor Q (X, m, r). The distribution of the 

efficiency factors at SAGE II aerosol wavelengths are displayed in Figure 

6.1. As can be seen, Figure 6.1 indicates that the best information of 

aerosol size distribution in the SAGE II aerosol extinction measurements is 

within the radius range between approximately 0.28 and 0.80 pm. A second 

source of the limits of the retrievable size range of a set of optical 

measurements is governed by the so-called ratio-criterion. According to this 

criterion, there is no information about aerosol size distribution that can 

be extracted from optical measurements at any two wavelengths in a radius 

range where the ratio of the two corresponding efficiency factors is a 

525 1 02 

constant. The calculated results of the ratio of Q * /Q * and 

e e 

385 ] 02 

/Q^ ‘ are given in Figure 6.2. Figure 6.2 also shows the ratio of 
.532 1.02 

Qg * * • As one can see from Figure 6 .2 , the nonretrievable radius 

range lies below the radius of 0.1 pm and above about 1.0 pm. In addition, as 
one would expect, the information content near these limits becomes 
relatively poor as revealed by the distributions of the efficiency factor in 
Figure 6.1. 
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eric aerosols 





6.1.2 Data comparison 


November 30, 1984 (Laramie, Wyoming) 

Figure 6.3 shows the stratospheric aerosol size distribution at 
altitudes (a) 15 km; (b) 18 km; (c) 21 km; and (d) 24 km. The two single 

mode aerosol size distributions are derived from SAGE II aerosol extinction 
measurements alone. The bi-mode model distribution is retrieved using both 
SAGE II and lidar aerosol optical measurements. Figure 6.3 indicates the 
close agreement among the three model distributions in the radius range 
between approximately 0.2 and 0.7 pm. Large variations are shown outside 
this range. This feature is consistent with the conclusions mentioned in 
Sect ion 6.1.1. 

Figure 6.4 shows the comparison of the backscatter coefficients obtained 
from lidar observations and from calculations based on aerosol size 
distributions retrieved from SAGE II aerosol extinctions. In general, there 
is a reasonable agreement between the observations and the derived results 
based on retrieved model size distributions, except near the altitude of 24 
km. 

A comparison of N and N 25 between dustsonde observations and 
calculations based on aerosol size distributions retrieved from SAGE II and 
correlative lidar observations is provided in Figure 6.5. Figure 6.5 
indicates that the retrieved model size distributions yield values of N 

• ID 

which are greater than those of dustsonde observations for all three analytic 
models. On the other hand, substantial improvement is found in the results 
of N 25 calculated using the bi-lognormal representation. 
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Figure 6.3b. Stratospheric aerosol size distributions derived from 
SAGE II aerosol extinctions at 18 km altitude. 
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Figure 6.3c. Stratospheric aerosol size distributions derived from 
SAGE II aerosol extinctions at 21 km altitude. 
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Figure 6.4a. Comparison of the backscatter coefficients obtained 
from lidar observations and from calculations based 
on aerosol size distributions retrieved from SAGE II 
aerosol extinctions. 

6-10 



ALTITUDE, KM 



BACKSCATTER COEFF, 1/(KM*SR) 


Figure 6.4b. Comparison of the backscatter coefficients obtained 
from lidar observations and from calculations based 
on aerosol size distributions retrieved from SAGE II 
aerosol extinctions. 
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January 9, 1985 (Laramie, Wyoming) 


Figure 6 . 6 a to 6 . 6 c are the same as Figure 6.3, except that they are for 
data obtained at altitudes of 15.2, 18.3, and 21 km, and that the dustsonde 
measurements of N N j^, and N 25 were used together with SAGE II aerosol 
extinctions to derive a bi-modal size distribution. In Figure 6 . 6 , the dots 
show wire-impactor data obtained by V. Overbeck of NASA Ames Research Center. 
Again, Figure 6.6 indicates the general agreement among the three different 
model size distributions in the radius range between approximately 0.3 pm and 
I.O pm, especially at an altitude of 15.2 km. In addition, at an altitude of 

15.2 km, the model size distributions agree well with the results of the wire 
impactor, except at r = .05 pm. 

6.2 Lidar Measurements in the Upper Stratosphere and Mesosphere 

The NASA-Langley 48” lidar system has been used for routine 
ground-correlative measurements in association with both the SAGE I and SAGE 
II satellites. The long term data set on stratospheric aerosol backscatter 
obtained with this instrument, which extends back to 1974, is a unique source 
of information on the long-term climatology of the aerosol in the northern 
hemisphere. Current use of the system is mainly confined to the region of 
the atmosphere below 30 km in altitude. This concentration of effort owes 
much to the fact that the stratospheric aerosol layer lies almost entirely 
within this region. The 48” lidar system was nevertheless designed with the 
inherent capability of being used to probe much higher in the stratosphere. 
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Figure 6.6a. 


Stratospheric aerosol size distributions derived 
from SAGE II aerosol extinctions at 15.2 km 
altitude. 
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Figure 6.6c. Stratospheric aerosol size distributions derived from 
SAGE II aerosol extinctions at 21.0 km altitude. 





with a measurable backscattered signal potentially obtainable from the 
mesosphere and even the lower part of the thermosphere. Although aerosols, 
in the form of noctilucent clouds, are known to exist in these higher 
regions, the main interest is in the variations in the molecular air 
density. These variations occur on a seasonal basis but, more importantly, 
are associated with the passage of atmospheric gravity waves, as well as 
tidal oscillations and turbulence. 

As presently constituted, the 48” lidar system cannot achieve its 
maximum potential for high altitude measurements. It is desirable to 

incorporate a more powerful laser and also to convert the current analog 
measurements of backscattered light intensity to photon counting. In order 
to assist with the instrument upgrade, the science requirements for high 
altitude work have been examined and simulations have been made of potential 
system performance. As a basis for quantifying the science requirements, 
literature pertaining to atmospheric structure and changes in the height 
range 30-90 km has been reviewed. This information has been obtained by 
lidar (Kent et al., 1972; Kent & Keenliside, 1974; 1975;Chanin & Hauchecorne, 
1981; Hauchecorne and Chanin, 1982; Chanin and Hauchecorne, 1984), rocket-borne 
instruments (Philbrick et al., 1980, 1983), and from Shuttle re-entry (Price, 
1983). These techniques show large density excursions to be common, 
particularly in the upper part of the indicated height range. Based on this 
information, two levels of desirable system performance were defined for the 
lidar system. These are shown in Table 6.1. The first level is that which is 
necessary in order to demonstrate the lidar system capability and to resolve 
the largest atmospheric changes occurring within the mesosphere. The second 
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TABLE 6.1. Science Requirements for High Altitude Lidar 


A. Minimum to Demonstrate Capability and Potential 

Lidar system must be able to make density measurements to the 
following specifications: 


Height Range 30-70 km 
Vertical Resolution 2 km 
Integration Time ^ 1 hr 

Altitude Required 
(km) Accuracy 


30 1% 

40 1% 

50 2% 

60 5% 

70 10% 


B. Minimum for Useful Scientific Measurements 


Height Range 30-90 km 

Vertical Resolution 1.5 km 

Integration Time 1 15 min 

Altitude Required 

(km) Accuracy 


30 1% 

40 1% 

50 1% 

60 1% 

70 2% 

80 3% 

90 5% 
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level is that required for routine scientific measurements to improve our 
understanding of atmospheric dynamics at these altitudes. It should be noted 
that not only do the two tables differ in the measurement accuracy required, 
but also in the altitude resolution and, in particular, the integration 
time. Short integration times of the order of 15 minutes are required to 
resolve the higher frequency gravity wave oscillations. It is also desirable 
that measurements should be made throughout a 24-hour period but limitations 
to present lidar technology will almost certainly restrict the highest 
altitude measurements to night-time only. 

Simulations of potential lidar performance have been carried out using 
the characteristics of the 48” lidar receiver, assuming photon counting, used 
together with one of two existing lasers or a third more powerful but 
potentially available laser. The results of these simulations are shown in 
Table 6.2 (using the existing Nd/YAG laser). Table 6.3 (using the existing 
Ruby laser), and Table 6.4 (using a more powerful Nd/YAG laser). Both Nd/YAG 
lasers are used at their second harmonics (0.53 pm) because of the much 
higher photomultiplier quantum efficiency available at that wavelength . 
Vertical resolution has been set at 2.0 km and integration times are either 
600 or 3600 seconds. The levels of system performance shown in Tables 6.2, 
6.3, and 6.4 may be compared to the science requirements shown in Table 6.1. 

If we first examine Tables 6.2 and 6.3, we may compare the relative 
merits of two existing lasers. Although the single shot performance of the 
Nd/YAG system is not as good as that of the Ruby lidar, the integrated 
performance is superior at all altitudes. The only caveat to be applied to 
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TABLE 6.2. High Altitude Lidar Simulation 


(A) SYSTEM PARAMETERS 

2 

Receiving Mirror Area - 1.00 m 

Full Angle Beam Divergence - 0. lOOE-02 r 

Wavelength - 0.6943E-06 m 

Overall Efficiency - 0.02 

Pulse Energy - 1.50 J 

Firing Rate - 0.13 Hz 

Filter Bandwidth -* O.lOOE-08 m 

Photomultiplier Noise - 50.0 

(B) ATMOSPHERIC PARAMETERS 

Atmospheric Transmission - 0.80 
Sky Brightness = 0.70E+00 R.A"1 

(C) OPERATIONAL PARAMETERS 

Height Increment - 1.0 km 
Integration Time - 3600.0 s 

System Constant - 0.224E+21 Sky Background = 0.117E-02 
PM Background - 0.667E-03 


HFIGHT 

SIGNAL 

NOISE 

4, UNCERTNTY 

r UNCFRTNTY 

(KM) 

COUNT 

COUNT 

(ONE SHOT) 

(INTEGRATED 

5 

.274E+07 

.183E-02 

• 060 

.003 

10 

.385E4-06 

.ie3E-02 

• 161 

.006 

15 

.806E+05 

.183P-02 

.352 

.017 

?0 

.150E + 05 

.1B3E-02 

.816 

.038 

25 

.507E404 

.183E-02 

1.294 

.061 

50 

.lQOF+04 

.183E-02 

2.292 

.108 

35 

.643E403 

.183E-02 

3.943 

.186 

40 

.232E+03 

.183E-02 

6.559 

.309 

45 

.904E+02 

.183F-02 

10.519 

.496 

50 

.3e2E+02 

.183F-02 

16.174 

.762 

55 

.175E402 

.ie3E-C? 

23.022 

1.126 

60 

.BCIE+Ol 

.1P3E-02 

35.346 

1.666 

65 

.350E401 

.183F-C2 

52.757 

2.487 

70 

.157E+01 

.183E-02 

79.777 

3.761 

75 

.660E400 

.183F-02 

123.217 

5.809 

PO 

.26BE+00 

.183E-02 

193.670 

9.130 

P5 

.1C6E+00 

.183F-02 

309,973 

14.612 


OO 

.393E-01 

.183E-02 

515.907 

24.320 

05 

.145E-01 

.183E-02 

882.017 

41.579 

00 

.532E-02 

.183E-C2 

1589.526 

74.931 
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TABLE 6.3. High Altitude Lidar Simulation 


(A) SYSTEM PARAMETERS 

2 

Receiving Mirror Area - 1.00 m 

Full Angle Beam Divergence - O.lOOE-02 r 

Wavelength - 0.5300E-06 m 

Overall Efficiency - 0-08 

Pulse Energy - 0.15 J 

Firing Rate - 2.50 Hz 

Filter Bandwidth - 0.100E-08m 

Photomultiplier Noise - 50.0 

(B) ATMOSPHERIC PARAMETERS 

Atmospheric Transmission - 0.70 
Sky Brightness - 0.70E+00 R.A“^ 

(C) OPERATIONAL PARAMETERS 


Height Increment - 2.0 km 
Integration Time - 3600.0 s 


System Constant - 0.314E+20 Sky Background - 0.467E-02 


PM Background - 0.667E-03 


IGHT 

SIGNAL 

NOISE 

r UNCERTNTV 

t UNCE 

KM) 

COUNT 

COUNT 

(ONE SHOT) 

( INTEGI 

K 

.113F407 

.53AF-02 

• 09A 

.001 

10 

.159E+06 

.53AF-02 

.251 

.003 

15 

.333F+05 

.53AE-02 

.5AP 

• 006 

20 

.620E+0A 

.53AE-02 

1.270 

.013 

25 

.2A7F+0A 

•53AE-02 

2.01A 

.021 

30 

.7P6E403 

.53AF-02 

3.566 

.036 

35 

•266E+03 

.53AE-02 

6.135 

.065 

AO 

.960E+02 

.53AE-02 

10.206 

.106 

A5 

.373E+02 

•53AF-C2 

16.367 

.173 

50 

.15PE+02 

.53AF-02 

25.169 

.265 

55 

.722E+01 

•53AE-C2 

37.233 

.392 

AO 

.331F+01 

.53AF-02 

55.032 

.580 

A5 

.lAPE+01 

.53AF-02 

62.210 

• 667 

70 

.650F+00 

.53AF-02 

12A.561 

1.313 

75 

.273F+00 

•53AE-02 

193.310 

2.036 

PO 

.lllE+00 

.53AF-02 

307. AA7 

3.241 

«5 

•A37E-01 

.534E-02 

506. A93 

5.339 

<)0 

.162E-01 

.53AE-02 

90A.351 

9.533 

95 

.598E-02 

.53AE-02 

1776.390 
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TABLE 6. 4. High Altitude Lidar Simulation 


(A) SYSTEM PARAMETERS 

2 

Receiving Mirror Area - 1.00 m 

Full Angle Beam Divergence - 0.100E-'02 r 

Wavelength - 0.5300E-06ni 

Overall Efficiency - 0.08 

Pulse Energy - 0.40 J 

Firing Rate - 10.00 Hz 

Filter Bandwidth - 0.100E-08m 

Photomultiplier Noise - 50.0 s“l 

(B) ATMOSPHERIC PARAMETERS 

Atmospheric Transmission - 0.70 
Sky Brightness - 0.70E+00 R.A’l 

(C) OPERATIONAL PARAMETERS 

Height Increment - 2.0 km 
Integration Time - 600.0 s 

System Constant - 0.837E+20 Sky Background - 0.467E-02 
PM Background - 0.667E-03 


MFT6HT 

SIGNAL 

NOISE 

Z UNCFPTNTY 

Z UNCERTNTY 

(KM) 

COUNT 

COUNT 

(ONE SWOT) 

(INTEGRATED 

*> 

.302F407 

.53AF-02 

.058 

,001 

10 

.A2AF+06 

.53AF-02 

.15A 

,002 

15 

,MR8F*C5 

.53AF-02 

.336 

,006 

?0 

.165F+05 

.53AF-D2 

,778 

,010 

25 

,658E-»0A 

.53AE-02 

1 .233 

,016 

30 

.210E-»0A 

.53AE-C2 

2,18A 

,028 

35 

.708F+03 

.53AE-02 

3.757 

.069 

AO 

.256E+03 

.53AF-02 

6,250 

,081 

A5 

.096E4C2 

.53AE-02 

10,022 

.129 

50 

.A21E+02 

.53AE-02 

15,611 

.199 

55 

.193E402 

.53AE-02 

22,795 

.296 

60 

.8P2E+01 

.53AE-02 

33,683 

.635 

65 

.3O6F401 

.53AE-02 

50,287 

,669 

70 

.173E401 

.53AE-02 

76,083 

.902 

75 

*■^28^400 

» 53A F-02 

117,666 

1.519 

«0 

•296E400 

.53AE-02 

185,551 

2,395 

P5 

.117E400 

•53AF-02 

299,636 

3.866 


90 

.633E-01 

.536E-02 

509.207 

6.576 

95 

.160E-01 

.536F-02 

916,612 

11.808 

00 

,5P6F-02 

.536E-02 

1805,101 

23.306 
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these results is in connection with the atmospheric transmission and, in 
particular, the background sky brightness. The values assumed correspond to 
good, low-light seeing conditions. Local measurements of these quantities at 
the 48” lidar site do not exist and may well not approach these optimum 
values. The performance figures for the Nd/YAG system may also be compared 
to the science requirement figures given in Table 6.1, It can be seen that 
level A (minimum to demonstrate capability and potential) is easily attained 
with a considerable safety factor. Level B (minimum for useful scientific 
measurements) is not attained (allow for different integration times), 
indicating the desirability of replacing the existing Nd/YAG laser by a more 
powerful model. Table 6.4, based on the performance of a better Nd/YAG laser 
may be compared with the figures for required accuracy in the second part of 
Table 6,1. The two tables have slightly different spatial and temporal 
integration widths (1,5 km versus 2 km and 600 S versus 900 S), the overall 
effect of this is, however, fairly minor. It can be seen that an upgraded 
48” system will meet the science requirements at all altitudes up to 80 km 
and comes close to the requirement at 90 km altitude (6.6% versus 5.0% 
accuracy). Once again, the caveat concerning atmospheric transmission and 
sky brightness must be considered. It should also be noted that these 

science requirements do not address the subject of daytime measurements. 
These are desirable scientifically and are possible up to a somewhat lower 
altitude given substantial modifications to the receiving optics of the 48” 
lidar system. Work carried out under a separate contract (NAS 1-17072) shows 
that on upgraded systems, using an interference filter to reject unwanted sky 
background, could make measurements with 10% accuracy (for 10-minute 
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integration period) up to about 60 km. Use of a narrow-band Fabry-Perot 
etalon might enable this limit to be extended up to about 75 km. 

6.3 Forty-eight Inch Lidar Receiver Modification 

As part of the work of modifying the 48” lidar receiver for the high 
altitude studies discussed in Section 6.2, a detailed analysis was made of 
the lidar receiver optics. The basic purpose of the analysis was to find a 
suitable way to incorporate a photon counting photomultiplier together with a 
high-speed rotating shutter. The latter is required to reduce, or eliminate, 
the signal induced noise in the photomultiplier caused by the initial, high 
intensity, backscattered signal from the lower atmosphere. The shutter, 
which screens the photomultiplier during the first 100-200 microseconds after 
the laser is fired is driven by a high-speed synchronous motor. Figure 6.7 
shows part of the existing optical system of the laser receiver, together 
with the proposed shutter/photon counting assembly. Collimating lens A is 
part of the existing system, lens B is required to focus the beam to a small 
diameter at the shutter itself. Figure 6.8 shows a scale drawing of a 

possible design for the shutter assembly. Final dimensions may change 
slightly, contingent upon the choice of shutter motor and any other changes 
made to the optical system. 

In the course of examining the optical system, it was noticed that the 
collimating lens (lens A in Fig. 6.7) was not of sufficient diameter to pass 
all nonaxial rays. Its diameter is 1.75”, exactly the same as that of the 
telescope exit pupil. Loss of light at this lens occurs (shown by the two 
outlying rays in Fig. 6.7) becoming increasingly serious for greater beam 
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Shutter Blade 
5.0" Diameter 



Figure 6.8. Diagram of a possible arrangement for the high speed shut ter /photon 
counting assembly. 


divergencies. That this loss actually took place was confirmed by 
measurement of the beam diameter on the telescope itself using the sun as a 
source and varying the size of the aperture stop- Correction of this fault 
appears at first sight to be simply a matter of replacing the collimating 

lens by one of equal focal length and larger diameter. This, although 

theoretically correct, raises other practical problems. The first is that 
the turning mirror would also have to be replaced by one of greater 
diameter. The second and more important problem is that, after passing the 
collimating lens, the beam passes through one or more interference filters. 
The standard diameter for these is 2.0”. For angular divergences greater than 
about 0.5 mr, these filters will also vignette the beam. A simple lens of 
shorter focal length placed closer to the field stop would solve these 
problems but the optimum location is at or very close to the turning mirror 
making the solution practically difficult. The solution chosen is to replace 
the collimating lens by two lenses as shown in Fig. 6.9. In this figure, the 

two new lenses and the path of an edge beam is shown as an overlay on the 

original diagram (Fig. 6.7). In the new optical arrangement, the system will 
pass all the light out to a beam divergence (total angle) of 2.0 mr. The 

collimating lenses have diameters of 2.0” and focal lengths of 20.0” and 
27.5”, respectively. The turning mirror does not need to be replaced and the 
second collimating lens is placed at the location of the current collimating 
lens. The resultant beam following the collimater is sufficiently narrow to 
pass through the interference filters and the entrance lens of the proposed 
shutter assembly. 
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Diagram showing possible modification to the 48” lidar 
receiver to correct for loss of radiation. Two new 
lenses are shown in black, the path of an edge ray by 
cross-hatch ing . 




7. TASK 6--INVESTIGATE STRATOSPHERIC PLANETARY HAVES AND THEIR EFFECT ON 


TRANSPORT OF AEROSOLS, 0^, HATER VAPOR AND NO^ 

7.1 Aerosols As Dynamic Tracers and Association with Planetary Have Activity 

Previous work carried out by Kent et al. (1985), has shown how the SAM 
II 1 pm extinction ratio may be used as a tracer for atmospheric movements 
associated with the north polar vortex. A similar study has been made on the 
use of the corresponding extinction ratios obtained from SAGE I and 
preliminary SAGE II data. Examination of extinction ratio profiles from 
different latitudes show that it is, in general, possible to classify these 
as a low-latitude type (layer maximum at about 25 km or above) or a middle or 
high latitude type (layer maximum at 15-20 km). At certain times of the 
year, however, some low latitude-type profiles are observed at high 
latitudes. Examples of these profiles obtained with the SAGE I satellite are 
shown in Figure 7.1. Figure 7.1(a) is an example of a "normal" middle or high 
latitude extinction ratio profile. Figure 7. Kb) shows a quite different 
profile, although it was obtained at the same latitude and on the same day. 
The latter profile is, in fact, similar to profiles normally seen within 
about 20° of the equator. Examination of 11 months of SAGE I data, taken 
prior to major volcanic eruptions, shows a very clear seasonal cycle in the 
occurrence of "low-latitude type" layers outside the equatorial belt. Table 
7.1 shows how the frequency of occurrence varies in both hemispheres with 
season. The occurrence is clearly a winter-time phenomenon; in summer, in 
both hemispheres, very few cases are observed. The most obvious explanation 
for this happening is that we are observing transport of low latitude aerosol 
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TABLE 7.1. Frequency of Observation of "Low Latitude Type" Extinction Ratio 
Profiles: February - December 1979, SAGE I Data. 
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into the winter hemisphere; this is a time of planetary wave activity and of 
the existence of the polar night vortex with its associated subsidence. 

The same phenomenon manifests itself in a similar form in the SAGE II 
data. The basic SAGE II extinction ratio profiles show much higher values of 
the extinction ratio due to the presence of material from the El Chichon 

eruptions. The "normal" layer at high latitudes has a peak extinction ratio 

of about 10 and a fairly flat layer top between 20 and 25 km. An example of 

this is shown in Figure 7.2(a). As in the case of SAGE I, anomolous layers 

also appear during the winter season. A typical example is shown in Figure 

7.2(b) where a secondary peak occurs at an altitude of about 27 km. These 

layers are observed in both hemisphere and, in as far as it is possible to 
see in the limited preliminary SAGE II data set, have a seasonal cycle 
similar to that observed with SAGE I. A detailed study has been made of the 

preliminary SAGE II extinction data between October 1984 and May 1985 and the 

following general characteristics of these anomalous layers determined. Part 
of this work has been presented at the Sixth Conference on Atmospheric 
Radiation, Williamsburg, Virginia May 12-16, 1986 (Kent, 1986). 

1 • Altiitude of Occurrence 

The effects are most visible at altitudes of 24-27 km. This is 
determined mostly by the relative shapes of the high latitude aerosol layer , 
where the "anomalies" appear , and the low latitude aerosol layer from which 

they are supposed to be derived. The exact altitude range is still 

uncertain , but there appears to be a change in behavior at an altitude of 

around 22 km. In order to illustrate this change. Figure 7.3 shows mean 
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Figure 7.2. High latitude SAGE II extinction ratio profiles, April 20, 1985. 


Potential, Temp. 



Figure 7.3. SAGE II extinction ratio profiles, November 4 - 
December 6, 1984. 
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aerosol extinction ratio profiles for two separate latitude bands during the 

period November 4 - December 6, 1984* The extinction ratio profiles are 

plotted in terms of potential temperature rather than altitude (the altitudes 

are marked on the figure). It can be seen that at all altitudes above about 

16 km, the low latitude aerosol has a much higher mean extinction ratio. The 

figure also shows a third curve derived from the spread in aerosol 

o o 

characteristics at 40 -50 . The meaning of this third curve is best explained 
by reference to Figure 7.4. This figure shows the probability distribution of 
extinction ratio values at an altitude of 25 km, the same latitude band of 
40^-50^N and for the period November 9-20, 1984. Most extinction ratios have 
a value of 4 or less, but there is a "tail” of high values extending up to an 
extinction ratio of 16. These values are the ones attributable to material 
transported from low latitudes. In order to describe this tail, we have 
chosen the 95 percent probability level as a representative value. This is 
shown in Figure 7.4 and also in Figure 7.3, as a function of potential 
temperature (and altitude). The interesting feature in Figure 7.3 is the 
slope discontinuity between 22 and 25 km. Two possible explanations exist 
for this change. 

a. Less material is transported from low to high latitudes at 
altitudes below 22 km as compared to the higher altitudes, or 

b. Material is transported at all levels but is more rapidly mixed 
into the background below about 22 km. 

2. Longitudinal Characteristics 

At or near 25 km, the occurrence of high extinction ratios is periodic 
with longitude. Wave 1 and Wave 2 structure is clearly evident and this is 
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EXTINCTION RATIO 


Figure 7.4. Probability distribution of extinction ratios at 
an altitude of 25 km. 
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correlated with the structure in the 30 mb pressure-altitude field. An 
example of the correlation is shown in Figure 7.5(a) and (b). Figure 7.5(a) 
shows the location of SAGE II measurements on November 13-15, 1984. The 

location of the largest extinction ratios is shown in black, the smallest by 
clear open circles. Figure 7.5(b) shows the 30 mb meteorological analysis 
map for November 14, 1984. The dominant features are the polar low and an 

elongated high pressure region over Japan. Comparison of 7.5(a) with 7.5(b) 
shows that the anomalously high values of extinction ratio occur to the north 
of the high pressure region where one might expect to find them on the theory 
of transport from low latitudes. This correlation is also evident over a 
longer period, as is shown in Figure 7.6. This figure shows the mean 
cross-correlation between the 24 km extinction ratio and the 30 mb pressure 
altitude for a period of about one month. The latitude range is approximately 
30° to 50°N. Maximum correlation of about 0.5 occurs at zero longitudinal 
shift. The periodic nature of the function is evident as is also a bias to a 
westerly longitudinal shift. Figure 7.6 shows the correlation function 
transformed to give the cross-spectral power density; almost all the power is 
in the first three planetary wave numbers. For the first two, there is a 
phase shift to the west. This is consistent with the occurrence of a net 
eddy transport of material from low to high latitudes by these waves. 

3. The Magnitude of the Extinction Ratios 

At the 25 km level between October and December 1984, the normal 
extinction ratio between 40° and 50°N was 4 or less. Anomalous values up to 
19 were observed. Such values are commonly observed at the same altitude in 
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Figure 7.5a. SAGE II extinction ratios, 24 km altitude, November 13-15, 
1984. 

Extinction O 8-12 

Ratio Key ^ 4-8 ^ >12 
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Figure 7.5b. 30 mb meteorological analysis map, November 14, 1984 
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Figure 7.6(a). Mean cross-correlation, November 7 - December 3, 1984 (24 km extinction 
ratio vs. 30 mb pressure altitude). 
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Figure 7.6(b). Cross-spectral power density, November 7 - December 3, 1984 
(24 km extinction ratio vs. 30 mb pressure altitude). 


low latitudes ( I0-20^N) during the same time period. (See Figures 7.3 and 
7.4.) 

4. Slope of Isentropic Surfaces 

Assuming that material is transported from low latitudes by planetary 
wave activity, it is likely that it will remain on a surface of contrast 
potential temperature. Thus, material origination at 20 km at 15^N will 
reach 45^N at an altitude of about 19 km, material at 25 km will suffer 
little height change, while material originating at 30 km will ascend to an 
altitude of about 31 km. The detailed relationship has already been shown in 
Figure 7.3. 

5. Seasonal Characteristics 

In the northern hemisphere, the phenomenon is limited to the period 
October 1984 - March 1985, the aerosol layer returning to normal in April, 
May, 1985. This agrees with the SAGE I data. In the southern hemisphere, the 
seasonal variation is similar but the period of absence of anomalous layers 
is much shorter, extending only from January through March, 1985. Again, this 
is in agreement with the SAGE I data and reflects the longer period of 
planetary wave activity in the southern hemisphere. It may be noted that 
SAGE I data indicate that there can be considerable year-to-year variation in 
both hemispheres. 

6. Relationship to Planetary Wave Activity 

A study has been made of the amount of aerosol transport as a function 
of time in the northern hemisphere from October 1984 - May 1985. Figure 7.7 
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shows histograms of the extinction ratio values for the latitude band 40~50^N 
for each satellite crossing period. Anomalous values, corresponding to 
extinction ratios of 4.0 or greater, are shown in black. The build up in 
activity from October to early January is clear; this was followed by a 
completely clear period, then a minor reoccurrence and a final disappearance 
in April 1985. The clear period in January is unexpected and corresponds to 
the time period following an unusual major sudden warming. A preliminary 
analysis of the variation in eddy transport during this period, made by 
Labitzke and co-workers at the Free University of Berlin (Labitzke et al., 
1985), shows a disappearance of planetary wave activity in the second part of 
January followed by a weak reestablishment of the polar vortex in 
February-March, 1985. This behavior, which is shown in Figure 7.8, 

corresponds very well with that of the satellite extinction data. 

7. Other Effects 

The recently discovered height error in the preliminary SAGE II sunrise 
data has injected an altitude uncertainty of the order of 1 km into that half 

of the analyzed data set. Re-analysis of this portion will be required after 

this error is corrected. 

Apart from the above, which is an error, there appears to have been an 
altitude change taking place over much of the northern hemisphere in the 

height of the top-side of the aerosol layers. This results in a steady drop, 
from October 1984 to January 1985, of the aerosol concentration at altitudes 
near 30 km. This effect requires further study. 
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Probability distribution of 1 pm extinction ratios at an altitude of 25 i 
and 40° — 50°N, October 1984 — April 1985. Anomalous values (extinction 
ratios >4) are shown in black. 
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Figure 7.8. Meridional time sections from November 1984 to March 1985 
of zonally averaged eddy momentum transport (m^/s^) and 
eddy heat transport (Km/s) by planetary wave 1. 

(Labitzke et aL , 1985.) 
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8. TASK 7--INVESTIGATE THE CONNECTION BETWEEN THE VARIATION OF AEROSOL 

BACKSCATTER AND LOCAL MICROPHTSICAL CONDITIONS 

8.1 Analysis of Data from the Airborne Lidar Flights of January 1984 

An important feature of the global stratospheric aerosol climatology is 
the occurrence of polar stratospheric clouds (PSCs) at high latitudes in 
winter. These were initially detected in data from the SAM II satellite 
(McCormick et al., 1982) and are also visible on occasion in SAGE I and SAGE 
II data. A limitation of the satellite extinction measurements is the 
restriction of the measurement location to the sunrise and sunset 
boundaries. It is thus not possible to study conditions within the polar 
night region where temperatures are colder and theoretical analysis (Steele 
and Hamill, 1981), leads us to expect a greater likelihood of PSC 
occurrence. In January 1984, the NASA-LaRC airborne lidar system made a 

special mission to high northern latitudes. The purpose of this mission was 


a. to make correlative measurements with the SAM II satellite, 

b. to map the volcanic aerosol injected by the eruption of El Chichon 
in April 1982, and 

c. to observe polar stratospheric clouds should arctic temperatures be 
low enough to permit their formation. 

Following a series of correlative measurements on January 23, 1984, 

three flights were made to the north of Thule, Greenland, including one to 
the north pole. Polar stratospheric clouds were seen on all three occasions. 
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on one occasion extending from bout 80 N to the pole. As part of the general 
analysis of stratospheric aerosol properties carried out in conjunction with 
the analysis of SAGE II data, a detailed study has been made of the lidar 
data from these flights. This work has been documented in the form of a 
research paper which has recently been accepted by the Journal of Atmospheric 
Sciences for publication. A copy of the abstract of this paper is attached 
to this report as Appendix 3 (Kent et al., 1986). 
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9. CONCLUSIONS 


A detailed study has been made of the preliminary SAGE II satellite data 
set, together with available associated data sets. The SAGE II data set 
extends from October, 1984 to May, 1985. Correlative data obtained during two 
measurement series made respectively on November 30, 1984 in Laramie, 
Wyoming, and on August 7, 1985 at Fairbanks, Alaska, have been analyzed in 
detail. In addition, climatological data on the stratospheric aerosol have 
been available from the SAGE I and SAM II satellite as well as ground-based 
and airborne lidar. 

In general, the SAGE II data set has been found to be both internally 
self-consistent and in good agreement with the results of other experiments. 
It is clearly a data set of great potential value. During the course of the 
study, various anomalies and inconsistencies in the data set have been found 
and removed from successive editions. Some still remain, although it is 
anticipated that these will be removed from the next version of the data set 
which should have archival status. Among those errors still remaining in the 
current preliminary data set are (1) an altitude error in the sunrise data, 

(2) anomalous values at low altitudes in the 0.525 pm aerosol channel, and 

(3) contamination from aerosol in the water vapor channel. The latter effect 
is clearly decreasing as the residual aerosol from El Chichon is removed from 
the stratosphere. The sources of these anomalies have been identified in the 
data inversion procedure and their correction is currently under study. 

In the course of examination of the SAGE II data, various climatological 
studies have been made in which the SAGE II data have been examined in the 
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context of other data sets. The agreement here is excellent. Detailed 
studies have also been made on the retrieval of the aerosol size distribution 
from the SAGE II multiwavelength aerosol data and its use in conjunction with 
the SAGE II water vapor channel to study the aerosol microphysics. 

Although the current study has been made on a preliminary data set and 
from that aspect should not be regarded as final, the general data quality 
appears extremely good. Known anomalies will be corrected in the next data 
version which will then be suitable for general scientific use and, as such, 
of unique value to many atmospheric studies. 
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ABSTRACT 


Satellite occultation measurements (SAGE I, SAGE II, and SAM II), 
and the results from recent airborne lidar expeditions, have led to a 
greatly improved understanding of stratospheric aerosol climatology. 
Marked seasonal and latitudinal variations are observed, particularly 
in the polar winter where polar stratospheric clouds occur and the 
polar vortex modifies the aerosol distribution. Several volcanic 
eruptions since 1979 have produced major changes in the aerosol con- 
centrations. 
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ABSTRACT 


Several eruptions of the El Chlcbon volcano between March 28 and April 4, 
1982 produced a global enhancement of the stratospheric aerosol content which 
seriously affected satellite measurements of properties of the earth and Its 
atmosphere. This enhancement, possibly the largest in the past seventy years, 
has since been monitored by in situ, as well as, satellite and ground-based 
remote sensing techniques. The injected material, at latitude 17 ”N and at 
altitudes up to 30 km, spread rapidly zonally and more slowly meridionally. 
Airborne Ildar measurements showed that the maximum global mass loading of 
about 1.2 X 10^ tonnes occurred 3-6 months after the eruptions. By May 1983, 
the aerosol was well distributed globally, with peak concentrations within a 
rather broad equatorial band and between 40^ and 60°S and S. The altitude of 
the peak aerosol mixing ratio, initially at about 27 km In the equatorial zone, 
had decreased by 5-6 km by May 1983. Hldlatitude Ildar measurements showed a 
similar decrease, accompanied by loss of material from the stratosphere Into 
the troposphere. Estimates of the 1/e decay time for the aerosol mass loading 
vary between S and 14 months, depending upon latitude and measurement 
techniques, with a northern hemisphere average of about 10 months. 
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ABSTRACT 


Airborne lidar measurements of backscattering at 0.6943 urn from polar 
stratospheric clouds, made in January 1984, are reported. The clouds, whose 
altitudes and geographical locations coincided with ambient atmospheric 
temperatures below about 193 K, were observed to cover a greater area of the 
polar cap than had previously been apparent from satellite measurements. 
They were seen on three separate flights north of Thule, Greenland (76.5°N, 
68.7^W), on one occasion extending continuously from approximately 80^N to 
the North Pole. Pronounced layering of the clouds was observed and the 
maximum backseat ter enhancement, relative to that from the background 
aerosol, was between 100 and 200. These values occurred at an altitude of 
about 20 km, close to the region of minimum stratospheric temperature. 
Depolarization of the order of 20-50% in the backscattered signal was 
measured, in support of the hypothesis that the aerosols forming the clouds 
are frozen. Comparison of the experimentally determined 

backseat ter ing- temperature relationship with a theoretical model , based on a 
VO lean ic aerosol and using best available estimates for water vapor 
concentration, shows good agreement at the 100 mb and 70 mb pressure levels. 
A sma 1 1 systemat ic error at the 50 mb and 30 mb levels may be due to 
inaccurate character izat ion of the temperature field at these altitudes and 
locations . 
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